Bose-Einstein Correlations in a Fluid Dynamical Scenario for Proton-Proton Scattering at 7 TeV 
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Using a fluid dynamical scenario for pp scattering at 7 TeV, we compute correlation functions for 
TT^TT^ pairs. Femtoscopic radii are extracted based on three-dimensional parametrizations of the 
correlation functions. We study the radii as a function of the transverse momenta of the pairs, for 
different multiplicity classes, corresponding to recent experimental results from ALICE. We find the 
same decrease of the radii with kr, more and more pronounced with increasing multiplicity, but ab- 
sent for the lowest multiplicities. In the model we understand this as transition from string expansion 
(low multiplicity) towards a three-dimensional hydrodynamical expansion (high multiplicity) . 
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I. INTRODUCTION 

Bose-Einstein correlations have proven to be a very use- 
ful tool to provide space-time information about colliding 
systems at relativistic energies. Different kinds of reac- 
tions have been considered, elementary ones like electron- 
positron annihilation [1], or more complex systems like 
proton-proton and nucleus- nucleus collisions [2. Sophis- 
ticated methods have been developed in particular for 
heavy ion collisions at 200 GeV, to interpret the depen- 
dence of the two particle correlation function C on the 
relative pair momentum q. A simple way to summarize 
the results amounts to parametrize the correlation func- 
tion as 
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where "long" refers to the beam direction, "out" is paral- 
lel to the projection of the pair momentum P perpendic- 
ular to the beam, and "side" is the direction orthogonal 
to "long" and "out" The fit parameters i?out, Rside, 

and i?iong will be referred to as "femtoscopic radii" in the 
following. In this way one can study for example the 
dependence of these radii on the centrality of the reac- 
tion, which is interesting because for central collisions we 
expect higher energy densities and finally more collective 
flow compared to peripheral collisions. A very general 
feature in heavy ion collisions seems to be the fact that 
all the radii increase with centrality. It is also found that 
for a given centrality all radii decrease substantially with 
increasing transverse momentum kx = |Pt/2|. This is 
compatible with a scenario of collective flow with a strong 
space-momentum correlation, as we are going to discuss 
later. 

Recently first results have been shown concerning the 
kx dependence of the femtoscopic radii, for different mul- 
tiplicity classes @. The amazing result: with increasing 
multiplicity, one observes a more and more visible de- 
crease of the radii with kx, as in heavy ions. And the 



radii increase with multiplicity, similar to the increase of 
the radii with centrality in heavy ion collisions. So do 
we have the same dynamics which governs proton-proton 
and heavy ion collisions? Is there a collective expansion 
driven by hydrodynamics in pp scattering as well? 

To contribute to the answer to these questions, we will 
postulate a scenario of hydrodynamic evolutions based 
on flux-tube initial conditions, and compute correlation 
function, make the same fitting procedures as in the ex- 
periment, and analyze the kx dependence of the femto- 
scopic radii. 



II. HYDRODYNAMIC SCENARIO 

Originally hydrodynamics was only thought to present 
a valid description for almost central collisions of heavy 
nuclei, where the volume is (relatively) big. But it seems 
that this approach works very well for all centralities. 
There is also no fundamental difference seen between 
CuCu and AuAu, although the copper system is much 
smaller. So it seems that systems much smaller than 
central AuAu fit well into this fluid picture. Finally it is 
more and more accepted that the famous ridge structure 
observed in angle-rapidity dihadron correlation [7] is due 
to fluctuating initial conditions, which are subsequently 
transformed into collective flow Here, the relevant 
scale for applying hydrodynamics is not the nuclear size, 
but the size of the fluctuations, which is typically 1-2 fm. 

Hydrodynamics is derived from a gradient expansion. 
In order to justify a hydrodynamical treatment one has 
to relate the length scale with the viscosity rj. There 
are many estimates of numerical values for the rj/S {S 
being the entropy density), however, all are based on un- 
proven model assumptions. As shown in fs^, the variation 
of the ratio of elliptical flow to eccentricity can be per- 
fectly explained based on ideal hydrodynamics, whereas 
other authors extract a non- vanishing viscosity from the 
same observable. It should also be mentioned that in 
the ideal hydrodynamical partonic phase is followed by 
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a highly viscous hadronic one, so in the average the sys- 
tem is viscous. By talking about ideal hydrodynamics, 
we mean an ideal hydrodynamical partonic phase, and as 
shown in [8*1, at least all RHIC AuAu data on soft physics 
are best described with viscosity zero. From a theoret- 
ical point of view, the viscosity of a QGP is unknown, 
and even a lower limit for rj/S as l/47r is not a mathe- 
matical limit, but rather an estimate with unknown error 
compared to the QCD value. 

Is QGP formation a nuclear phenomenon? Or can it 
be formed in pp scattering, as already proposed in 
[T^ . [isl [T^ ? Based on the above discussion, there is no 
reason not to treat proton-proton scattering in the same 
way as heavy ions, namely incorporating a hydrodynam- 
ical evolution. This approach makes clear predictions for 
many variables, so the Nature will tell us whether the ap- 
proach is justified or not. Therefore it will be extremely 
interesting to think about the implications of such a mini 
QGP, how such a small system can equilibrate so quickly, 
and so on. It would be an enormous waste of opportuni- 
ties, not to consider this possibility, since a vast amount 
of proton-proton data will be available very soon, con- 
cerning all kinds of observables. 

What makes pp scattering at LHC energies interesting 
in this respect, is the fact that at this high energy multiple 
scattering becomes very important, where a large number 
of scatterings amounts to a large multiplicity. In such 
cases, very large energy densities occur, even bigger than 
the values obtained in heavy ion collisions at RHIC - but 
in a smaller volume. Several authors discussed already 
the possibility of a hydrodynamical phase in pp collisions 
at the LHC, to explain the ridge correlation or to predict 
ehiptical flow [13 18]. 

We are going to employ a sophisticated hydrodynami- 
cal scenario, first presented in ref. |8| where many details 
can be found, with the following main features: 

• initial conditions obtained from a fiux tube ap- 
proach (EPOS), compatible with the string model 
used since many years for elementary collisions 
(electron-positron, proton proton), and the color 
glass condensate picture |2C)| : 

• event-by-event procedure, taking into the account 
the highly irregular space structure of single events, 
being experimentally visible via so-called ridge 
structures in two-particle correlations; 

• core-corona separation, considering the fact that 
only a part of the matter thermalizes |21|]; only in 
the core region, the energy density from the strings 
is considered for the hydrodynamical evolution; 

• use of an efficient code for solving the hydrodynamic 
equations in 3+1 dimensions, including the conser- 
vation of baryon number, strangeness, and electric 
charge; 
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Figure 1: (Color online) Expansion of a low multiplicity (a) 
and a high multiplicity (b) event. 

• employment of a realistic equation-of-state, com- 
patible with lattice gauge results - with a cross-over 
transition from the hadronic to the plasma phase 

msi; 

• use of a complete hadron resonance table, making 
our calculations compatible with the results from 
statistical models; 

• hadronic cascade procedure after hadronization 
from the thermal system at an early stage (zj . [25| . 

In ref. jSj] , we test the approach by investigating all 
soft observables of heavy ion physics, in case of AuAu 
scattering at 200 GeV. In refs. [l^, [l^ we investigate 
first proton-proton results, with among other things "the 
ridge". 

III. RESULTS 

In the following, we will consider several multiplicity 
classes, named mult 1, mult 4 , mult 7 and mult 8, corre- 
sponding to four out of the eight multiplicity classes used 
in ref. [6], going from low multiplicity {mult 1, less than 
minimum bias) to high multiplicity {mult 8, five time min- 
imum bias) . Our core-corona procedure will find no core 
for mult 1, then increasing core fraction, and for mult 4 
to mult 8 essentially core, with increasing energy densi- 
ties. So the mult 1 events are just ordinary strings which 
expand longitudinally (see fig. HJa)), whereas mult 4 to 
mult 8 show a hydrodynamical expansion, also in trans- 
verse direction (see fig. (Hb)). In fig. [51 we show the evo- 
lution of the energy density for the different multiplicity 
classes. Obviously the energy density increases with mul- 
tiplicity, values of more than 100 GeV/fm^ are achieved. 



Bose-Einstein Correlations in a Fluid Dynamical Scenario for Prototon-Proton Scattering at 7 TeV 



3 




1 2 1 2 1 2 

r (fm) r (fm) r (fm) 



Figure 2: (Color online) Evolution of the energy density (aver- 
age over many events) for the three multiphcity classes Milt 4, 
Milt 7 and Milt 8. The full red lines are the initial conditions, 
the other curves are respectively 1 and 2 fm/c later. 

We consider two options for the equation of state, one 
being a parametrization of the results of [l^l , the other 
one referring to [23'|. There is a big difference between 
the two, the transition temperature is much lower and 
the transition is much smoother in[25|, compared to [23]. 
The corresponding "freeze out" radii (where the transition 
hydro / cascade is done) are very different, see fig. [31 The 
final results will differ less, because the early freeze out 
for the EoS from [22] is followed by an intense hadronic 
rescattering. 

Based on roughly ten million simulations of the hy- 
drodynamical evolution, we compute in the usual way 
the correlation functions for 7r+7r+ pairs, taking into ac- 
count Bose-Einstein statistics, as discussed in (sl. [Tol [26|. 
Whereas in the data Pythias has to be used as "baseline", 
we stay consistently within our scenario and use simply 
a calculation without Bose-Einstein statistics as baseline. 
We then fit the correlation functions to obtain the radii. 
Before showing the results, let us discuss in a qualitative 
way why we expect a decrease of the radii with fcy. Let 
us consider the freeze out curves of fig. [3] Most particle 
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Figure 3: (Color online) Freeze out radii, using the EoS from 
[2a | (solid curves), and using [2^ (dashed curves). 



production occurs in the region where the radii drop to 
zero. Comparing the two points "1" and "2" in the figure, 
we have to recall that the collective fiow at a large radius 
(1) is much bigger compared to small radius (2). In fig. |4j 
we sketch the corresponding momentum vectors of pairs, 
emitted at large and small radii (which finally amounts to 
large and small kx), where the vectors are such that their 
difference is the same. The space distances in case 2 are 
then bigger than the ones in case 1 (and these distances 
are essentially the femtoscopic radii). 

Let us discuss our main result: As seen in fig. [SJ all 
radii indeed show a more and more pronounced decrease 
with increasing k^, for data and simulations, which can 
- in the calculations - clearly be attributed to collective 
flow. For the case Milt 1 the radii i?out and i?side are es- 
sentially flat, only i?iong has already some kx dependence. 
So we see here nicely the transition from a longitudinal 
expansion (string) towards a three-dimensional hydrody- 
namical expansion for higher multiplicities. 
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Figure 4: (Color online) Emitting pairs with fixed momentum 
difference in radial direction, from different radii. 
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Figure 5: (Color online) Fenitoscopic radii for three different multiplicity classes, using an EoS compatible with ref. (2^ (solid 
curves) or compatible with ref. [23| (dotted curves). The curves are the same for Milt 1 because there is no fluid phase. The 
points are data. The curves are absolute predictions, the parameters of the model are obtained from other comparisons (yields, 
pt spectra). 
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